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ABSTRACT 


This work is a finite element analysis of an external braking 
system , to determine contact stresses and slip . The study is 
restricted to a static case and the body forces are neglected . To 
limit the study to two elastic bodies m contact , lever is 
eliminated from the analysis and suitable force boundary 
conditions are applied to the top surface of the shoe . The 
standard two — body algorithm , which has a fast convergence rate 
is employed „ The results are compared with the relations used in 
the conventional brake design procedure . The effect on contact 
stresses and slip are studied for change in parameters like the 
total contact angle , the coefficient of friction , the hinge 
location etc . 



Chapter !_ 


Introduction 


1 . 1 Brakes 


Man has gained speed in all aspects of life . Be it a 
vehicle , machinery , or any other equipment , high velocity is a 
necessity . But they cannot run continuously as a perpetual motion 
machine . So a device to control their motion is needed . Brake is 
a device which by means of frictional resistance is used to either 
retard or stop a rotating member . 

Brakes are classified on the basis of actuating mechanism 
as Hydraulic C Hydrodynamic Brakes and Fluid Agitator 3, Electric C 
Generator and Eddy Current Brakes 1 and Mechanical C Shoe , Band 
and Block - Band 3 . Mechanical brakes are further classified into 
two groups according to the direction of braking force 

(1) Radial Brakes C External and Internal 3 

(2) Axial Brakes C Disc and Cone 3 

i 

In this thesis, only mechanical braking system witl* 
external shoe is considered . A typical system is shown in fig l.i 
• The essential part is a lever , with fulcrum at "F" to whicl 
shoe “S" is bonded . By applying external force "W" at the end o 
the lever , the shoe is brought in contact with the drum " D “ an 
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b r ak i ng *r' f i or pi ^cp - » 

fig. 1,1 . When the drum is rotating in clockwise direction the 

free body diagram of the lever ~ shoe assembly is shown in fig- 
1.2 . "N" is the resultant normal force and " tM u is the frictional 

resistance due to braking . 

As the moments of frictional force and the external load 


,, W“ have the same direction , the frictional force is assisting 
the external load in the operation of braking . this type of brake 
is called self - energizing . If the direction of rotation in fig. 
1.1 is reversed , the frictional force force will oppose the 
externally applied load as far as the moment about fulcrum is 
concerned . This is a de — energizing system of braking . 

When the moment of all the forces about fulcrum is equated 


to zero , we get 


W » 


N <b + l/z-/ja ) 


( 1 . 1 ) 


For b >> 1/2 


W = 


N (b— £ia) 


( 1 . 2 ) 


When (b-/L/a) < 0 , M W 4< becomes negative . This situation is 
called self - locking . In self - locking brakes , external load 
is not applied to brake the drum . The shoe is merely allowed to 
come in contact with the drum . The frictional force generated by 


the contact itself is sufficient to brake the system . 


A brake should be self - energizing to reduce the external 






load , But it should not be self - locking to avoid stoppage of 
the drum if , by accident * the lever fails on the drum . 

The factors T which a f ~~ ~ c v t * * * p e ** f q r mane e of a braking 
system are as follows : - 

(1) Coefficient of friction between the braking surfaces. 

<2) Ability of shoe material to dissipate heat . 

(3) Area of contact surface . 

(4) Average and maximum pressure on the braking or contact surface 

(5) Rate of wear . 

(6) Velocity of the drum periphery . 

To control the above factors for the favorable outcome , 
the major governing criteria are the selection of the shoe 
material and the braking or contact area of the shoe . The 
material with a high dynamic coefficient of friction is normally 
selected so that it can generate the resisting torque necessary 
for the braking action . 

When a brake operates , heat is generated due to the 
presence of friction and slip at the drum — shoe interface . The 
generated heat depends on the coefficient of friction , the 
pressure distribution on the braking surface , the velocity of 
drum periphery , the inter facial slip and the area of brakinc 
surface • Since the temperature rise affects adversely thi 
strength , the coefficient of friction and other characters of th| 


shoe material , it should have high heat dissipating capacity 



Sometimes it becomes necessary to provide some cooling 
arr angemen t . 

The b r a k i n g or the contact area of me shoe is selected so 
that the maximum pressure does not exceed the safe value for the 
shoe material . The maximum pressure also depends on the rate of 
wear . Velocity of drum periphery is , in general , not under the 
control brake designer but is a given parameter • 

1 ,2 Brake Design Procedure 

The primary purpose of a brake is to absorb energy and 
then dissipate it without permitting any appreciable temperature 
rise that can damage either the shoe or the drum surface . The 
brake design procedure , similar to designing any other mechanical 
device , starts with the establishment of general requirement for 
which the brake is to be designed . The primary requirement is the 
generation of braking torque sufficient to balance the drum torque 
for the given velocity • 

The traditional method of brake design uses the following 
equation to relate the braking torque "T" with design variables 
like the coefficient of friction "fj" , the maximum pressure "ptnox 11 
, the radius of the drum "r" , the shoe width "b" and the arc of 


contact ( Juvinall 


C13 ) 



C I j 



Here B is the angle that a radial line on the drum makes 
with the line joining the drum center to the fulcrum and S± and Bz 
are it's values at the contact points < fig* 1.3 ) . This formula 

is based on many idealizations • It is assumed that the pressure 
on the braking surface is proportional to the normal wear 
Further , while calculating the normal wear , it is assumed that 
everything except the brake lining has infinite rigidity . This 

gives the following distribution of pressure along the contact 
surface . 


P 


pmcLX 


< sin0)max 


(1.4) 


First a suitable material for the shoe < and the drum if 
the selection of the drum material is m the jurisdiction of the 
brake designer ) is chosen . Once the material is selected , the 
coefficient of friction between the mating surfaces and the 
capacity of the shoe material in terms of maximum allowable 
pressure "pmax " are known from the data books . Then a suitable 
fulcrum for the brake lever is chosen .Finally for the given drum 
of radius "r" and width H b" , the arc of contact ( Bt - Bz > is 
calculated iteratively from eqn. 1.1 . 


The value of /j and p given in the data book are usually 
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for normal operating conditions 9 when it is assumed that the heat 
generated due to friction will be dissipated . Since both th* 
coefficient of friction as well as the capacity to withstand 
pressure decrease with temperature , it is important to estimate 
the temperature at the friction - surface during braking and 
calculate the value of and *‘p l£ correspond ing to the actual 
temperature • Unfortunately , it is quite difficult to calculate 
the temperature . 

When the designer doesn't have temperature data available 
, he must rely on the principle that the rate of wear is 
proportional to the rate at which energy is absorbed by the 
friction surfaces ( Phelan C23 > . This is because the temperature 
is also proportional to the rate of energy absorbed , which using 
the principle t the capacity of the lining material is expressed 
in terms of the product ppV y the rate of energy absorption per 
unit area . Dubey and Suresh C33 have suggested an empirical 
approach for two - wheeler's brake design in which the contact 
area and the drum diameter are related to the torque based on the 
study of the two - wheelers using new technology • 

It is quite clear that the assumptions on which the eqn. 
1.3 is based are too simplistic . For an optimum brake design , 
the rate of wear and the temperature distribution at the friction 
surface should be evaluated more accurately . Both these 
quantities depend on the rate of friction work which in turn 



depends on the pressure distribution and the slip at the friction 
interface - But: the pressure distribution and slip depends on the 
coefficient of friction and the contact geometry which are 
influenced by temperature and wear . Because of this the inter 
dependence , the procedure has to be iterative and incremental 
It is obvious that for a brand new brake for which the wear is 
zero , the pressure and slip at the friction surface will be 
governed solely by the material properties of the shoe and drum , 
the coefficient of friction and the initial contact geometry 
Once the pressure and the slip are known at the friction surface , 
the temperature at the initial instance can be determined using 
appropriate conduction equation . It can be used to modify the 
material properties and the coefficient of friction for the next 
iteration . The iteration can be continued till we get the 
convergence values of the pressure , slip and temperature « Next , 
the rate of wear at the initial instance can be calculated using 
appropriate diffusion equation . Using this information , we can 
modify the contact geometry for the calculation of the next 
instance • 

In this thesis , the problem of determination of pressure 
and slip only for the brand new brake is undertaken . It is 
treated as a contact problem between two elastic bodies with I 
friction • The external loading is both normal and tangential and 

i 

the initial contact surface is circular • 




1.5 t i t e rature Survey 

Several investigators have attempted to determine stresses 
and displacements in two elastic bodies having a curved contact 
surface . In three dimensions , the problem of contact of two 
elastic bodies was first formulated and solved by Hertz C43 under 
several restrictive assumptions • Research work based on Hertz's 
theory up to 1982 has been reviewed by Johnson C53 . Hertz assumed 

that the geometry of general curved surfaces in the vicinity of 
contact can be described by quadratic terms only . Then the 
geometric compatibility condition leads to an elliptical contact 
area . Neglecting friction and assuming that the bodies in contact 
deform as though they were elastic half ~ spaces , Hertz used the 
Boussinesq solution to deduce that the contact pressure 
distribution has to be elliptic to produce an elliptic contact 
area . Hertz evaluated the dimensions of the contact area and the 
stresses at the contact surface but presented only a speculative 
sketch of the sub surface stresses . The stresses beneath a generai 

elliptical contact were analyzed first by N. M. beiajet in two 

\ 

papers in 1917 and 1929 (see ref. 5) • It is observed that maximum^ 

shear stress occurs at a depth of (0.57a) for a circular contact: 
of radius ' a ' . 

From the Hertz's theory of contact stresses , where the! 
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contact is over a finite area f the case of two - dimensional 
contact between two cylinders can be obtained as a limiting case 
by allowing the contact area to become infinitely long in one 
direction . However , it is easier to do it directly as shown by 
Poritsky C63 . He also evaluated the internal stresses due to 
tangential traction at the surface . He observed that , for a 
static contact , slipping occurs outside the central zone of 
length (2b # ) which is given by 

b * = b ^1 - Q / (fjP) (1 .5) 

Here b is the semi - contact length , fj is the coefficient of 

friction and Q and P are respectively the tangential and normal 

loads such that Q < ( /_iP ) . On the other hand , for a rolling 

contact , the sticking region is no longer centrally located but 

is adjacent to the leading edge of the contact « 

For a three - dimensional sliding contact (Q = fjP ) f the 

internal stresses due to combined effect of normal pressure and 

tangential traction were obtained by Bryant & Keer and Sack field & 

Hills (see ref. 5) . For a static contact , Cottaneo (see ref .5) I 

showed that , to avoid slip the tangential stress at the boundary! 

of contact area has to be infinite . Thus , some slip is 

t 

i 

inevitable at the edge of contact area even under the action off 
the smallest tangential load . Like a two - dimensional problem , j 

A 'I 

here also the no — slip zone is a central region . For an, 

t 

elliptical contact area of semi — axis a and b , the semi - axes' 
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a"' ana D f of the no - siip ellipse are given by 

a' =aC 1-Q/ </LiP ) D 1X9 , b' =bC 1-Q/ <£/P ) 3 1X * (1.6) 

where tu , Q and P are as defined by eqn. 1.5 . The annular zone of 
'micro - slip' penetrates further into the contact area with 

increasing Q until , when Q=/jP the whole contact slides . Being a 

dissipative process , the micro - slip has the effect of making 

the state of contact stress dependent upon the history of 

i 

loading. j 

The effect of friction is to bring the point of maximum shear 

j 

stress closer to the surface . For values of fj ,0.3, the maximum 

| 

shear stress occurs at the surface . 

Another important analytical method to solve contact problems 
(with or without friction) is to formulate it in terms of integral 
equations . Muskhelishvili C73 and Gladwell E03 solved some 
contact problems using this approach . Other notable solutions, 
which use the method of integral equations are those of Weitsmari 
C93 for unbounded contact between a plate on the elastic half H 

space , Prasad & Dasgupta C10D for rectangular block compressed by 

] 

rigid planes and Spez*ce till for indentation of elastic half - 
space by an axisymmetric punch . [ 

When contact geometry cannot be described adequately by the 
quadratic terms f it becomes necessary to use numerical methods 
Two approaches have been proposed • In the first approach , the 
conditions of displacement compatibility are supplied to ari 
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assumed contact area and the equilibrium equations are solved to 
find the contact stresses * If near the edges ? the normal 
stresses come out to be tensile , that region is excluded from the 
contact area in subsequent iterations . Further , if' the shear 
stresses exceed the product of coefficient of friction (fu) and the 
normal stresses (tn> , the tangential compatibility condition is 
replaced by the slip condition in subsequent iterations . The 
iterations are continued till the normal stresses become 
compressive and the shear stresses become less than or equal to 
/lit n . The second approach is based on the principle that the 
pressure distribution and the area of contact will be such as to 
minimize the total elastic strain energy , subject to no 

| 

interference outside the contact and positive pressure inside . | 

t 

When friction is present at the contact , one needs an additional 
minimization principle to find the slip . The advantage of the 
method is that the process of minimization can be carried using j 

t 

mathematical programming . 

Many investigators have adopted the technique of finite 
element method to solve the elastic contact problems . Here only a 

5 

* f 

few typical papers are reviewed . The papers which use the first 

t 

approach are discussed first . Ohte C123 has proposed a method in | 

l 

f 

which the stiffness matrix is modified to eliminate the reaction 

forces at the contact nodes based on stick — slip contact j 

! 

conditions . On the other hand , Sachdeva and Ramakrishnan C133 f 
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havp used the flexibility matrix to satisfy the frictional contact 
conditions • The flexibility matrix is obtained by inverting the 
stiffness matrix , which is first condensed to eliminate all the 
nodes except those which are likely to be in contact and those 
where external forces act . Gautam et al C143 have compared the 
method of Ohte and Sachdeva & Ramakr ishnan and concluded that m 
problems involving no symmetry of geometry or loading , the second 
method is more economical in terms of memory and solution time 
requirements . Whereas proportional loading has been assumed in 
both these methods , Gaertner C 15 3 has used incremental loading 
Triangular elements with 6 dof (normal displacement , tangential 
displacement , rigid rotation , tangential strain , normal stress 
and tangential stress > have been used . Use of stress components 
as dof facilitates adoption of any frictional law and not just 
Coulomb's law . Elasto-plastic frictional contact problem has been 
treated by Valllappan et al C163 with incremental loading 
Schreppers et al C173 have presented a finite element formulation 
for a large sliding contact . For such a problem , criterion to 
identify sets of actual contacting material points constitutes an 
important aspect of the contact algorithm . Dandekar and Conant 
E183 have used the boundary element method to solve frictional 
contact problem. 

In the above — mentioned works , the contact conditions have 
been applied by modifying either the stiffness matrix or the 
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flexibility matrix - However , in the last dozen years o r ' so , ir. 

most of the papers , the contact conditions have been treated as 
additional constraints and lagrange multiplier or penalty function 
method have been used to satisfy them . Ok a mo to and akazawa E19D 
have used the method of lagrange multiplier to solve frictional 
contact problem with incremental loading . Chaudhary and Bathe 
C203 have developed a general algorithm for three - dimensional 
dynamic contact problem using lagrange multipliers . Stadter and 
weise C213 have developed a penalty method with the use of 'gap' 
elements where the gap element stiffness is taken as the penalty 
number . The stiffness is varied to account for regions of contact 
and separation . A gap element with friction was used by 
Mazurkiewicz and Ostachowicz C223 . Mottershead et al E233 have 
reviewed the literature on the use of lagrange multiplier or 
penalty function methods to frictional contact problem with 
incremental loading and presented a general finite element 
approach for the treatment of contact problems . The mathematical 
foundation of lagrange multiplier and penalty function methods by 


finite 

element technique 

is studied 

in ICikucHi and 

Odert 

C243 

Unlike 

the most problems 

in elastic 

stress analysis 

, the 

contact 

problems are posed as 

variational 

inequal i ti es 

owing 

to the 


unilateral form of the displacement constraints . 

The method of mathematical programming to non - Hertzian 
friction less contact problem was applied first by Kalker and 
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Hung and Sauxe C2&3 nave usea finite 


element method 


along with mathematical programming technique to solve friction 
less contact problem .Recently f Mahmoud et al C273 have used 
incremental convex programming (along with finite element 
technique to solve friction less contact problem . There does not 
seem to be any solution of frictional contact problem using the 
technique of mathematical programming . 

All the above papers have placed emphasis on presenting a 
general methodology with a few illustrative problems rather than 
solving engineering problems of practical interest . Needless to 
say that there are no studies of the contact stresses in brakes 
either analytical or numerical. 

1 . 4 Objective And Scope of Presen t Work 

In a brake shoe -drum assembly , cm$ contact surface is 
convex (i.e. that of drum ) while the other is concave ( i.e. that 
of shoe ) . Therefore one cannot use Hertz's solution . In fact , 
it would give infinite contact length and zero contact pressure as 
the radii of curvature are equal and opposite • As stated in the 
section 1.3 , there is no analytical solution . Therefore , a 
finite element method is used for det ermination of the contact 
pressure and slip at the shoe-drum interface . The problem is 
idealized as two — dimensional (plain stress) . The distribution 
of the external loading on the shoe is calculated from the 
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assumption that it is proportional to the deflection of the brake 
lever as it rotates to bring the shoe in contact with the drum 
The deflection is calculated from the geometrical compat lb i 1 i ty « 

The iterative algorithm proposed by sachdeva and 
Ramakr ishnan C133 has been used expect for one difference 
Whereas the contact condition have been applied to flexibility 
matrix in reference C133 , here they have been applied to 

stiffness matrix only • Like in reference C133 , the size of the 
global stiffness matrix has been reduced by skyline storage and 
condensation . The iteration algorithm converges quite fast and 
gives nodal displacements on convergence - The difference in the 
tangential displacements of the shoe and the drum at the contact 

! 

surface gives the slip . To find the contact stresses , first the 
stress tensors evaluated at the gauss points are extrapolated to 
find their values at the boundary points . Then the stress vector 
for each body is evaluated and the magnitude of the normal and !j 

tangential components of the two stress vectors are averaged . I 

I 

The contact pressure distribution and the maximum value 
are compared with the corresponding equations of the traditional 1 
method of brake design . The parametric study is carried out to 
study the effects of contact angle , the coefficient of friction 
and the location of hinge on the contact pressure and the slip . 

The present analysis is applicable to new brakes more j 

i 

effectively when wear is less . It can be extended further to find i 
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the wear analysis . Distribution cf external forces on the shoe 
boundary has been derived from the geometry of She - snoe 

assembly . The process nowhere uses any incremental loading 
procedure , where final result depends on the accuracy with which 
the intermediate results converge - Instead full external load is 
acting on the shoe from the first iteration . 

1 .5 Plan of Thesis 

Finite element formulation and sub structuring are 
discussed in chapter 2 . It also gives procedure for calculation 
of stress • Chapter 3 deals with d iscr l t isat ion , boundary 
condition and results with discussions . The conclusion and scope 
for further work are presented in chapter 4. 



( 20 ) 


CHAPTER 2 

Mathematical Formulation Of Contact Probl em 

In this chapter , finite element algorithm for elastic bodies 
in contact is presented . The main feature of the algorithm are 
sub structuring of the conventional stiffn-ess matrix to reduce 
it's size and the modification of the condensed stiffness matrix 
to take care of the contact conditions. Finally , the procedure 
for calculation of the contact stress is discussed as a post - 
processing phenomena . 

2 . 1 Finite Element Formul at ion 

A body can be specified by the domain Q surrounded by the 
boundary r . The whole boundary of the body is divided into two 
parts 

<a) Boundary where displacements are prescribed (Fu) 

(b)Boundary where surface traction T are prescribed (Tf) 

Fig. 2.1 shows a plate - like body of thickness h . For such a 
body plain stress situation prevails and stress - strain relation 
in vector form is given by 
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for plane stress 

E = Young's modulus of elasticity 
v = Poison's ratio 

Let the whole body is divided into four nod ed quadrilatera 
elements , in which node numbering scheme is as shown in figur 
2.2 « For a any element displacement at any point can b 
expressed in* terms of displacements at the nodes of the element 
i . e . 


< u >*= C N 3 < u > ne 


( 2 . 2 ) 


where { u >*= displacement vector at any point of a element 

r O O T 
= < U U > 
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and C N 3 = matrix containing shape functions 

= j^l ° ! N 2 ° ! ! N 4 0 1 


[0 N 1 : 0 N 2 ! ! o N 4 J 

Fig. 2.2 also shows the node numbering scheme for a particular 
element . For each element , local coordinate system (if, 77) is 
chosen such that the local coordinates of the nodes 1, 2, 3, and 4 
are respectively (-1,-1) , < 1 , — 1 > , (1,1) and (-1,1) . Then shape 
functions are defined as 


n i = -J- (1 -o<i-77> 

N 2 = ~4~ (1+? > (1_r J > 

N 3 * — <1+0 <l+7?> 

N 4 = ~~ <1-0 <l+7?> 


(2.3) 


The global coordinates x and y of a point in an element are 
related to the local coordinates and 7? by the equation 


< x y° « C N 3 < X > 


no 


(2.4) 


Where < X >° “coordinates of a point in an element 


= i x y > 


C X > n41 ^coordinates of the nodal point of the element 

T 

= ( x y x y x y x y > 
i 7 t 2 7 2 a -'a a * 

C N 3 = Shape function as defined in equation 2.2 



(2 


In finite element formulation „ relationship between strain 
displacement is expressed as : — > 


_ <5ux 

X = <5x 
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By using eqn. 2.2 it becomes 


< £ y = C B 3 C U > * 

r 6 


where 


C B 3 = 


<5x 

0 

6 


Sx 

S 


C N 3 


( 2 . 6 ) 


Sy Sx 

The potential energy of the system is the difference of tht 
strain energy U and the work W of external force . Thus 

n = U - W (2.7) 


where 


U - -±- .f , 
2 J n < 


£ y € cr y h dA 


and 


W = f < u > T C T > h dl 

*Vf 


Here the body forces have been neglected . Consider a typical 
element 'e' whose area is and which has the boundary Tf* asl 

i 

I 

a part of Tf . Note that for internal elements and for elements 

< « t 

which do not have any boundary as Tf has no contribution in W .! 



( 2 5 ) 


Substituting eqns • (2 - 1 ) and (2.6) into (2-8) « we get the total 

potential energy for the element ' e ' as 

= J Cu> n ®CB3 T CD3[B: Cu >"* dA 

Cia 

— T 

f < u > ne C N 3 T < T > dl (2.8) 

J (Tf ) 

For the element to be in equilibrium , the potential energy should 
be minimum - Setting the first derivative of rte to zero , we get 
C K 3° < u > ne = < F > e (2.9) 

where 

C K 1** = stiffness matrix for one element 

CB3 T CD3CB3dA 

L2g 

and i F = load matrix for one element 

= J rf [ N ] T { T } dl 

« 

In order to evaluate C K 3* by numerical integration , it needs to 
be transformed to coordinates by using the transformation 

eqn. (2.4) . Then dA becomes equal to | J | d? dr) , where | J | is 
the determinant of the Jacobian matrix given by 

[ <s*_ <*y i 

<5* <5t? 1 

<s* _gx_, J 

<5? 


JACOBI MATRIX 
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Here C K J ~ has been evaluated by 2 X *2 gauss - legendr** numerical 
integration scheme . After calculating the elemental stiffness 
matrix , it's elements are transferred to the global stiffness 
matrix , which is stored in the skyline form • To accommodate 
distributed loading on the surface , < T > on each side is 
expressed m terms of nodal traction values i T > n * . Thus 

{T>°=CN]{T> n * (2.10) 



The elemental load vector is also evaluated numerically by gauss — 
legendre integration scheme after the change of variable from 
(x,y> to (?,Tl> - After calculating the load vectors for all the 
elements , they are assembled into the global load vector . 


2.2 SUB STRUCTURING 


For contact problems , two bodies have to be analyzed 
simultaneously . This makes the size of stiffness matrix quite 
large , which requires more memory space on the computer as well 
as more computational time • To reduce the memory size as well as 



( 2 ?; 


computational time ? condensation of global stiffness matrix is 
carried out . As a first step m the processes of condensation the 
nodes are divided into three categories 

(1> Nodes where nonzero external load is applied or nodes which 
are in contact • They are called nodes of type 1 and are numbered 
first . 

(2) Nodes for which the specified displacements are zero . They 
are called nodes of type 3 and are not numbered . 

(3) All nodes except those lying in the first two categories are 
called nodes of type 2 and are numbered after 1st types of nodes . 
They include internal nodes as well as those external nodes where 
the specified external forces are zero . 

Consider the bodies 'A' and 'B' shown in fig* 2.3 . The 


application of the formulation of section 2.1 leads to 


C K A 3 C U A > = £F A > 


C K B 3 £ u®> = £F B > 


( 2 . 12 ) 


Partitioning C K A 3 according to the type of nodes ,discussed above 
we have 
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As <L u _ > -£05 , the first two set of equations can be written 


as 


rcK^cK^mcu^i rcF^>j 

[ C K^KK- 2 3j[<u^j = |<F->j 


<2. 14) 


Similar equations are obtained for body 'B' also . Combining these 
equations ,we have 
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(2.15) 


Considering the second equation and eliminating the braces for 
ease in writing, following equation is obtained : — > 

u 2* “ K 22 K 21 U 1 <2.16) 
putting the value of u^in equation number 1 of (2.15), following 

equation is obtained 


where 


K u i * F l 


C K 3 = C K -K K" 1 K 

11 12 22 21 


(2.17) 


3 



(2 9 ) 


Thus the final matrix size is reduced to that of 1st types of 

nodes . The displacements of type 2 nodes can be calculated using 

eqn. (2.16) whenever needed . Note that while evaluating C K 3 

from eqn. (2.17) ,one has to extract C K 3 , C K 3 and 

C K 1 = C K 3 T from the matrices C K A 3 and CK B 3 which are in 
12 2 1 

skyline form . 

2 . 3 Con tac t Cond i t ions And Their Appl icat ion 


All the nodes except those in the possible contact region 
have either the known force or the known displacement , But at the 
contact nodes , neither force nor displacement is known . So the 

A 

sub structured stiffness matrix C K 3 , obtained in the previous 
section , cannot be inverted to get the displacements . To 
overcome this difficulty the sub structured stiffness matrix is 

I 

modified with the help of contact conditions . As each contact | 
node on the first body has a corresponding node on the second body 
and each node has two degrees of freedom , there are four 
equations for every node on the contact surface . So to make the 
stiffness matrix solvable , four contact conditions are needed . 

The contact conditions are as follows t— 

(1) The normal forces at the mating contact nodes are equal and 
opposite . 

(2) The tangential forces are equal and opposite at the in - 



C 3 0) 


contact nodes - 

(3) If the nodes in contact are 

(a) sticking, the tangential displacements are equal and in the 
same direction . 

(b) slipping ,the tangential forces at the node should be fj 
times the normal force , Where is coefficient of friction . 

(4) The difference between the normal displacements is equal to 
the clearance between the contacting surfaces . 

To apply the contact condition to the bodies shown in fig. 
2.3 , where node 'i' of body 'A' comes in contact with the node 'j' 
of body 'B' . <2i)th and <2i-l)th rows of C K A 3 and <2j>th and 

A B 

(2j-l)th rows of C K 3 have to be modified . The details are 
described below . 

(1) The mathematical statement of the first contact condition is 
Fni = Fnj ,where nt and nj are the normal directions 
at the nodes ' i ' and 'j' . In terms of (x,y> components the 

condition becomes 

—Fix Sinot + Fiy- Cosa — Fjx Sinot + Fjy- Coset = 0 
where 'a' is the angle made by the tangential direction si with 

A. 

the x— axis . To apply the condition (2i-l)th row of K is replaced 
by 

_(2i-i)th row X Sinot + (2i)th row X Cosot 

- (2 j-1 > th row X Sinot 
and F (2 j— 1 ) is made zero . 


(2 j ) th row X Cosot 



(2) The second contact condition can be stated as - Fjs =0 

or in terms of the (x,y) components 

—Fix Cosot — Fly Sinot 4 Fjx Since 4 Fjy Si not = 0 
To apply this condition <2i)th row of C K 3 is modified by 
<2i— l)th row times (-Cose*) 4 (2i)th row times (—Smct) 

4 (2 j — 1 ) th row times < Cos a) 4 (2j)th row times < Sina) 
and F(2j) is made zero • 

It is to be kept in mind that after applying first contact 

A 

condition , <2i-l ) th row of C K Dgets modified • So (2i -1 ) th row 

A 

of C K 3 has to be stored for the purpose of applying the second 
contact condition . 

(3) The third contact condition has two forms . 

(a) If the node is sticking uls 4 ujs = 0 ,or 

uix Cos a 4 uiy Since - ujx Cosoe - ujy Since =0 . This condition is 

A 

used to replace (2j“-l)th row of C K 3 . Thus all the elements of 
this row are made zero except the following ( 

(2i— l)th element is made (Cosoe) , (2i)th element is made (Since) , 
(2j“-l)th element is made (-Cosa) and (2j)th element is made 
(-Sina) . Finally F(2j-1> is made zero . 

1 

(b) If the node is slipping , jFs| - /j |Fn| = 0 For the | 

l 

second body , this condition becomes I 

Fjx (ss XCosct -*■ sn X X Sinct) + Fjy <s« X Sinot - sn X (J X Coset) -O j 

f 

.where ' ss' and 'sn' are respectively the signs of Fjs and Fjn 

A 

To apply this condition [ K i is replaced by 





(ss X Cosa + sn X /i X Si net) times (2 j — I > th row * 

| 

(ss X Sina - sn X /u X Cosa) times (2 j ) th row . The right hand side 
is made zero . 

(4) The fourth contact condition uln - ujn = cUaranc« or 
“■uix S i not + uly Cosot — ujx S i not + u jy Cosot = cl«<xr<xr*c« is used to 

a 

replace (2j)th row of C K 3 . Thus all elements of the row are 
made zero except the following 

(2i—l)th element is made (-sina), <2i)th element is replaced by 
(Cosot) , (2j~l)th element changes to (-S met) and (2j)th element 
becomes (Cosot) . The right hand side is made equal to the 
clearance . 

A 

After modifying the stiffness matrix C K 3 and the load 
matrix f F > as described above , the displacements of the nodes 
of type 1 can be obtained by solving the modified system . However 
, the number of contact nodes i.e. the contact length is not known 
a priori • Therefore the set of equations (2.17) have to be solved 
iteratively as described in the next section • 

2.4 Solution Procedure 

1 A certain number of nodes are assumed to be in contact 
Further all of them are assumed to be sticking • The contact 
conditions are described in section 2.3 are applied to modify the 
matrix t K 3 and C F > • The modified system is solved to find the 



displacements of the nodes of type 1 . Then the nooa I forces at 
these nodes are calculates . 

2 The nodal forces are checked to determine the status of 
nodes for the next iteration . 

If the normal force at a node becomes tensile ( outward from the 
surface ) ,it is deleted from the possible contact zone for the 
next iteration . 

For the nodes remaining m contact the value of iFsj^jFnJ is 
compared with (j . If jFsjy^jFnj is greater than (J , it means the 
node slips and in the next iteration , the slipping contact 
condition C 3(b) of section 2.3 1 is applied at these nodes . 

3. The iterations are continued until all the normal forces 
in the contact zone became compressive and the ratio of tangential 
force to the normal force at all the contact nodes became equal to 
or less than the coefficient of friction . 

Once the iteration are over , the stresses are computed as 
described in the next section . 


2.5 Contact Stress Calculation 

It is observed that when bilinear approximation is used for 
displacements , stresses are accurate at the center of the element 
. Therefore , first the stress components are determined at the 


center of elements in the layers adjacent to the contact surface . 



The calculation of stress components is done using the eon. €7 J ^ 
and eqn - (2-6) - Then the stress components at boundary points are 

obtained by extrapolation . For example , to find the stresses at f 
point 'P' . ext rapolat ion is done from the center of the element 
'E±' and 'Ez' (see fig. 2.4) . Then the surface tractions at the f 

boundary points are determined from the equation 
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(2.18) 


where x> and x> are the direction cosines of the outward normal 
x y ' 


boundary 

. Then 

the normal and 

tangential stress 
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The same procedure is adopted to find the normal and tangential 
stress component on the boundary of the second body . Finally 
magnitude of the normal and tangential stress components are 
ajpra* ged over the two bodies • 
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CHAPTER 3 


RESULTS AND DISCUSSION 


In this chapter , analysis of an external shoe brake system 
is carried out as a plain stress contact problem using the 
formulation of chapter 2 . The variations of normal and shear 
stresses as well as the relative slip are obtained along tne 
contact surface . Comparison with the expression used in the 
conventional design procedure as well as parametric studies are 
also undertaken . 


5 , 1 Model inq Of The Problem 

The brake system is idealized as a contact problem between 
the shoe and the drum only . Since the lever has been completely 
eliminated from the analysis an important aspect of the modeling 
is the estimation of the distribution of the forces exerted by the 
lever on the shoe . To determine this distribution , we assume the 
following : — 

(i> The lever is rigid compared to the shoe . 

(ii) As the lever rotates about the fulcrum to make the contact 
with the drum , the normal stress exerted by the lever on the drum 
is proportional to the vertical deflection of the shoe - 



Consider the Geometry of the lever shown in the f iq. 3.1 . Ulher 

the rigid lever rotates through 6a about 'F' , the deflection oi 

the point 'Q ' on the 1 ever is p { 6cx) , where 'p' is the 
distance of the point from 'F' . If ' cx ' is the angular location oi 
the point with the horizontal , then the vertical components oi 
the deflection becomes p(<3cOCosc* . Since pCosot = b + x , the 
horizontal distance of 'Q' from the fulcrum , the vertical 
deflection 6v of the point becomes proportional to it's horizontal 
distance from the fulcrum , i.e. 

<5v °c ( b + x ) (4.1) 

Since the normal stress exerted by the lever is assumed to b«? 
proportional to 6v , it also becomes proportional to ( b + x ) 

Thus j the variations of the normal stress becomes linear . Th* 
modeling of the actual variation by a linear function seems to b« 
a better approximation than by a constant function ( i.e. uniform 

distribution ). To find the actual variation , one has to include 

. f j 

the lever and consider it as a multiple — contact problem . Thitj 

1 1 

is likely to make the problem very complicated . It is felt that k 

|| 

for a thick shoe , the derivation of the actual variation from thfy 

ii 

linear case may not affect the far away stresses at the shoe \ 

i 

drum contact , very much . Therefore , it is decided to use thf; 
linear variation for the normal stress . The distribution of th| 

I 

tangential force is assumed to be uniform . The bette| 
approximation for this distribution may be obtained by relating if 




not 


dons 


to the horizontal deflection of the shoe , but it is 
again felt that it may not affect the contact stresses very much 
If ft and f 2 are the values of the normal stress at the tw 
end points as shown m fig. 3.1 , then for the case of lines 
variation 


(4.2) 


Further , if 'N' is the value of the normal force pr unit sho, 
width required to generate the specified braking torque , then 


1/2 (fl+f2)l 


(4.3) 


Solution of the eqns. (4.2) and (4.3) gives 


fi = 


2 N b 

(2 b + 1 ) 1 


f 2 = fl + 


(4.4) 


2 b + 1 


The distri.^ tion scheme for the fulcrum and the shoe is show- 

in the fig. 3.2 . To solve the finite element equation of th 

# i 

i 

assembly , all the rigid body modes have to be constrained . Fd, 

| ! 

the drum , it is done by assuming that the part , which is ij ; 

contact with the shaft has zero displacements . Since , a fixeij 

;,i 

support gives rise to reactive forces , constraining any othev 

tj 

point will not be realistic . Further , this part of the drum i| 
far away from the contact with the shoe . So it will not affecJ 
the contact stresses appreciably . The compatibility of the normal 






d i sp I ac t bet sen the drum and shoe rs^ovss sc^5 of the r 2 d 1 d 

body modes of the shoe « However , the tangential force on the 
shoe is equal to fj times the normal force • Therefore 9 the shoe 
is force to rotate about the drum axis • To constrain this 
rotation , the slip at the right - most contact point of the drum 
i.e. the actual slip at a point minus the slip of the right most 
point . To find the absolute slip , we need to add the slip of the 
right most point , which can be calculated from the velocity of 
the drum periphery. 


5 « 2 Mater 1 al Proper 1 1 es And Other Data 

The following data , which is typical to a brake system is 
used in getting most of the results • In parametric studies the 
contact angle ( <p ) ? the coefficient of friction (^) and the hinge 

location (b) are changed . 

Material properties 

Drum cast iron . 

Young modulus of elasticity Ei = 1.0 X 10 N/mm 
Poisson's ratio v ■ * 0. 25 

Shoe • — sintered material 

Young modulus of elasticity E 2 = 6.0 X 10 9 N/mm* 

Poisson's ratio v - 0.3 

z 

For the sintered material on cast iron 




the dynamic coefficient 



of friction vanes from 0.15 to 0.45 , whereas the capacity for 

maximum pressure varies from i.030 iM/mm* to 2.070 N/mm 2 . The 

following values are used . 

Coefficient of friction ^ = 0.3 
Normal force on the brake N = 100.0 N/ntra 
Dimensions 

Drum radius r = 125.0 mm 
Shoe - thickness = 30.0 mm 
Shaft diameter = 10.0 mm 

Total contact angle of shoe and drum , <p - 60° 

Lever dimensions Cfig. 1.11 
a = 50.0 mm 
b = 200.0 mm 
c = 250.0 mm 

The drum is discritised using 168 elements and the shoe has 72 

elements . 13 nodes are initially assumed to be in contact . 


3 . 3 Contact Stresses And Slip 

Fig. 3.3 shows the normal stress variation over the contact 
surface . As expected f the variation has a peak to the right of 
the mid point of the shoe . To compare the variation and the value 
of the maximum normal stress with the corresponding relations used 
in the commercial brake design , eqn 1.3 is plotted in fig. 3.4 
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depends 


Oi 


the 


lever 


a ! ono with 


Lilts . As the 


&L|U6 it 1 US 5 


dimension a (see fig. 1.1) , it is plotted for several values of 

a in a certain range . As stated in chapter 1 , this equation 

has been derived by assuming the pressure to be proportional to 


the wear and by neglecting the elasticity of the a 


rum m 


estimating the wear variation . On the other hand our results are 
for a new brake for which the wear is zero and which take into 
account the elasticity of the drum . Further , in our formulation 
the pressure is independent of the lever dimension 'a' . So it is 
expected that these two graphs should not match . However , two 
points need to be made . First is that according to eqn < 1 f 
' p 4 doesn't drop to zero at the last contact point . This is not 
realistic . Secondly for typical values of 'a' , eqn. d.4> 
underestimates the value of pmox for a new brake . Since this 
value of pmox is used in choosing the other brake parameters in 
the conventional method , the method seems to give very 
conservative design . 

The variation of shear stress over the contact surface is 


shown in fig. 3.5 . Since all the contact nodes slip , the shear 
stress should be fj times the normal stress . Thus the variation of 
the shear stress should be the same as that of the normal stress . 
By comparing the fig. 3.3 and 3.5 , it is found that the variation 
is more or less same . But at some points , the values of shear 
stress exceed the product of {u times the normal stress by a 



-e.180 
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considerable amount . This ^ t 

<nis fli »Ci cpancy may be attributed tc 

numerical error as well as to fh. *,_4. ■ . , , 

wen as to the fact that the coulomb relation 

was applied to the nodal forces rather than the stresses . 

The variation of relative slip xs plotted in fig. 3.6 . It is 

observed that the relative slip is maximum at the left - most 

contact point and decreases gradually to zero at the right most 

contact point . Fig. 3.3 and 3.6 together show that , except at 

the end - points , product of the pressure and the slip is more or 

less uniform . Thus , at least for a new brake , the rate of wear 

will be uniform over most of the contact region . 


3.4 Parametric Studies 
( a ) Contact Angle 


Fig. 3.7 shows the variations of the normal stress alone 
the central surface for various contact angles . The pattern ai 
the variation doesn't change with the contact angle but the valut 

f' 

of maximum normal stress decreases with the contact angle . This. 

V 

variation is shown in fig. 3.8 . This is expected , as smaller 
contact angle means smaller contact area . But since the total L 
normal force remains the same , the maximum value has to increase 
. Fig. 3.9 shows that although the variation of slip doesn ‘ 

change with the contact angle , it's value increases 
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In 


c on v s n t i on 3 1 ces 


ion procedure 


once 


the material 


selected , the contact area is chosen from the condition that ptmxx 
should not exceed the capacity of the material to withstand 
certain pressure . From the above results , we can conclude that 
it is disadvantageous to choose a contact length longer than the 
required , as it would increase the frictional work considerably . 


<b ) Coef f ici en t of f rict i 


on 


Varying the coefficient of friction has very little effect 
on either the variation or the values of the normal stress . The 

1 

shear stress values , however , increase proport lonal ly with \j 

1 

Further , the values or the variation of slip also do not change 

appreciably * 

If it is possible to increase the coefficient of friction j 
of the shoe material without changing it's capacity for sustaining 
a certain pressure , then it is desirable to do so . In that case 
one needs smaller contact area to generate the same braking torque 

r 

m 

(c? Hinge Location 

Fig. 3-10 shows the variation of the normal stress along 
the contact length for three different hinge locations . The j 
parameter used to characterize the hinge location is 
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horizontal distance of the left edge of shoe from the fulcrum 
The first graph is for one extreme case of p=0 (i.e. triangular 

variation of normal stress) while the third graph is for the other 
extreme case of b— >inf inity (i.e. uniform variation of the normal 
stress) . The variation of slip for these three values of ' b ' is 
shown in fig. 3.11 . 

While the variation of slip is not affected by change in 
'b' , that of the normal stress does . For a triangular loading 

(b=0> , the variation of the normal stress is also triangular . It 

is observed that for b>>200 , there is hardly any decrease m the 
value of maximum normal stress . So no advantage is gained by 
increasing b beyond this . Instead , the space requirement of the 
braking system would be more . 



CHAPTER 


Qi=LD c I yjILLQn And Suggest i on For Future Work 
4 , 1 Cone lusion 

The following conclusions can be drawn regarding the 
analysis of brake - drum assembly as a two body contact problem 
(1) The algorithm which so far has been restricted to only 
straight contact surface and normal loading can be applied to 
curved contact surface and tangential loading as well . 

< 2 ) The slip at the contact surface is not uniform . However the 
product of the normal stress and slip is more or less uniform over 
most of the contact surface . Therefore , the wear will be uniform 

(3) Increasing the whole contact angle increases the amount of 

slip * 

(4) The normal stresses and slip are unaffected by change in the 
coefficient of friction . 

(5) Changing the distance of the fulcrum from the drum center 
beyond a certain value has no effect either on the normal stresses 
or slip . 



x i on hor ruture u/orx 1 

The following suggestions are made for future work 

( a) The analysis can be improved by including the body - forces 
due to rotation of the drum 

(b) The application of variable coefficient of friction , which ' 
depends on factors like the rise in temperature of shoe , the drum 

, * i 

velocity , the wear of shoe etc. ,will give a better estimate of $ 
maximum contact stresses or slip . 

1 

(c) Since the exact stress boundary condition on the top part of j 

the shoe are difficult to estimate , it is better to include the ' 

! 

lever in the analysis to make it a three — body contact problem - j 
To analyze such a problem , the present contact algorithm will j 
have to be modified . j 

<d) In our formulation , the normal stresses can depend on the I 
wear only through the change in contact geometry of the shoe . To 
take care of this dependence , one can develop an incremental 
algorithm in which the wear is calculated at periodic intervals of 
time and used to update the contact geometry of the shoe * 
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